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Abstract: Using a generalized coordinate along with a proper invertiblc coordinate transforma- 
tion, we show that the Euler-Lagrange equation used by Bagchi et al. [la] is in clear violation of 
the Hamilton's principle. We also show that Newton's equation of motion they have used is not in a 
form that satisfies the dynamics of position-dependent mass (PDM) settings.. The equivalence be- 
tween Euler-Lagrange's and Newton's equations is now proved and documented through the proper 
invertible coordinate transformation and the introduction of a new PDM byproducted reaction-type 
force. The total mechanical energy for the PDM is shown to be conservative (i.e., dE/dt = 0, unlike 

J_h , Bagchi et al.'s [16J] observation). 
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on: 

The growing interests and/or research developments in the position-dependent mass (PDM) quantum mechanical 

systems described (mainly) by the von Roos Hamiltonian [l| (see also the sample of references |2h12| and related 

references cited therein) have inspired the relatively recent and rapid research attention in the PDM for classical 

O |' mechanical systems (cf, e.g., [12H14I . Ilq . [La . KL7I) . In fact, the interest in the PDM classical particles dates back to 

' 1974 through the Mathews and Lakshmanan [18| study of the equation of motion 

C3 ; (l + £x 2 )x-£xx 2 +uj 2 o x = 0; £eR, (1) 



o 

CO 
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where an overhead dote indicates a time derivative. This equation corresponds to a nonlinear oscillator, exhibiting 
amplitude-dependent simple harmonic oscillations [19j, with the Lagrangian 



j?: l -2{-tt^) (2) 

1 1 and the canonical linear momentum 

p = m(x)x; m (x) = (3) 

1 + £ar 

In their attempt to reproduce the Lagrange equation in (1), Bagchi et al. fnf have proposed that in the absence of 
C*") . any external force term the Newton's equation of motion with PDM gets modified to 

m (x) x + m (x) x = 0, (4) 

o 

where the prime indicates spatial derivative. Which when compared with (1) would suggest a PDM function in the 
form of 

rn(x) = — , (5) 

b ■ 

on ignoring, of course, the presence of the harmonic term io 2 x to effect such a comparison. They have also recollected 
Cruz et al's work [14| and used 

d ( dL\ dL ~ ~ 1 / . . ., 

d-t{-0x)-lTx= R -> R =-2 m ^ X > (6) 

where R is a reaction thrust (as named by Cruz et al. (HI)- This would make their Lagrange's and Newton's equations 
of motion equivalent and consistent. However, this is an improper approach for the PDM-settings in the classical 
mechanical framework. 
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A priory, the Lagrange's equation, used by Cruz et al 14|, in (6) is a clear violation of the textbook Hamilton's 



least action principle, one of the most profound results in physics. Yet, one would notice that the Newton's equation 
in (4) is based on the conservation of the linear momentum, dp/dt = p = 0, which is an improper approach for PDM 
settings. Mazharimousavi and Mustafa [17] have very recently shown that the quasi-linear momentum (as maned 

therein) 



II (x, x) = y/m (x)x (7) 

is the conserved quantity (i.e., II (x, x) = IIo (xq, xq) and II (x,x) — , where xq and xq are the initial position and 
initial velocity of the PDM, respectively) and not the linear momentum (i.e., p (x, x) ^ po (lo, xq), and p (x, x) ^ 0). 
In this communication, we fix this issue through the following arguments. 

Let us consider a classical particle with "unit mass" in the generalized coordinate q moving with velocity q in a 
forcc-frcc field, V (q) = 0. Then the corresponding Lagrangian L (q, q) is given by 

L = L{ qi q) = \q\ (8) 



and satisfies the Euler-Lagrange's equation 



±(BL\_BL 

dt\dq) dq [ > 



Under such assumptions, Eq.(9) would that q = 0. Which is in exact accord with Newton's law of motion for a free 
classical particle moving in the generalized coordinate q with a conserved generalized momentum 

dq 
Pg =- = 0=^q = qo , 

where % is the generalized initial momentum. Now, let the coordinate transformation 



q = q(x)= / VfW)du => q' (x) = y/f(x) =^q(x) = x^f(x) (10) 



represent a mapping from the coordinate q onto the coordinate x. As long as our choice of the coordinates is 
invertible (i.e., det {dxi/dqi) ^ 0, in general, which is the case under consideration here) then one can easily show 
that the Euler-Lagrange equation in (9) reads 



d (dL\ dL 
dt\dx) dx 



dx\ d fdL\ dL . . 



This in fact follows immediately from the Hamilton's least action principle. The Euler-Lagrange equations are known 
to be coordinate invariant for they assume the same form in all coordinate systems, provided that the coordinates 
invertibility is secured. As such, the Lagrangian L (q, q) in (8) would transform into 

L(x,x)= l -f{x)x 2 (12) 

with the corresponding Euler-Lagrange equation 

f(x)x+^f'(x)x 2 =0. (13) 

One may now go backwards and start with the Langrangian L (x, x) — ^m (x) x 2 of a PDM particle m (x) = / (x) 
moving in a force-free field, V (x) = 0, and use (11) to obtain 

1 



m (x) x — ——m (x) x =^> \J m (x)x = \/m (xo)xq => H (x, x) — Hq {xq, Xq) . (14) 

This result would not only support our thesis in |17| on the conservation of the quasi-linear momentum (i.e., II (a;, x) = 
0) but also suggests the amendment that has to made for Newton's equation of motion (4) for PDM used by Bagchi 
ct al. [16] and by Cruz et al. [14| in the absence of any external force term. In what follows we shall see that the 
"any external force" terminology holds true only for constant mass settings. 



In this regard, we may express the linear momentum in terms of the quasi-linear momentum 



p = m (x) x = y/rn {x)H (x, x) , (15) 

and cast Newton's equation of motion as 

F ex t (x, x) = F ext (x) + Rp DM (x, x) = — , (16) 

at 

where F ext (x) = —dV (x) /dx = is the set of potential energy driven external forces, and Rpdm (x, x) represents 
any feasible PDM-byproducted reaction- type force (of course, if it exists at all). To find out such a PDM-byproducted 
reaction force, Rpdm (x,x), one would naturally find 

T~ = T t m (• E ) x] = m (x) x + to (x) x , (17) 

Li b IX h 



and 



dp d " 



dt dt 



1 m (x) . . . 1 



Rpdm (x, x) = -7- = — y/m (x)II (x, x) = -— ^===ill (x, x) = -m (x) i 2 . (18) 



2 \Jm (x) 2 



When (17) and (18) are substituted in the Newton's equation of motion in (16), along with F ext (x) = 0, one obtains 

—m (x) x = m (x) x + m (x) x =$■ m (x) x 4- — rn (x) x =0 (19) 

which is in exact accord with the Euler-Lagrange result in (14). Obviously, the equivalence between the Euler- 
Lagrange's and Newton's equations is now documented through the above "good" invertible coordinate transformation 
(i.e., dx/dq ^0^ dq/dx) and the introduction of the new PDM-byproducted reaction-type force Rpdm (x,x) into 
Newton's law of motion. Hereby, we may safely conclude that the PDM setting is nothings but a manifestation of 
some "good" invertible coordinate transformation that leaves the corresponding Euler-Lagrange equation invariant. 
Under such textbook documentation settings, one would write the PDM Lagrangian as 

L = T-V=-m(x)x 2 -V(x) , F ext (x) = -ii jt 0, (20) 

2 ox 



with the corresponding Euler-Lagrange equation 



!_'^ -2 , dV(x) 



m (x) x + -to (x) x z H -^ = 0. (21) 

The Mathews and Lakshmanan [18J equation of motion in (1) is obtained if a PDM particle with to (x) = (l + £x 2 ) 
is subjected to move in an oscillator potential V (x) = ^to (x) w 2 x 2 or a hypothetical potential V (x) = —\m (x) u; 2 ,/^ 
|20j (if such potential exists at all). Both potentials yield 

to (x) x — to (x) £xx + to (x) uj x = 0. (22) 

Therefore, in the presence of an external periodic force with additional damping term and a quartic potential, the 
result in (22) would immediately suggest the following amendment to Eq. (11) of Bagchi et al.'s [16j to read 

m (x) x — m (x) £xx +to(x) ut q x + Ax + ax = / cos uit. (23) 

Which for a constant unit mass (i.e., £ = 0) recovers a forced, damped Duffing oscillator. Yet, with x = y Eq.(12) of 
Bagchi ct al. [16] should (taking z = ui) read 

. = fox 2 -^ + (1 + ^ 2) [fmsz _ Xx3 ay] 
1 + £x z 

Moreover, one can show that the time rate of change of the total mechanical energy E = T + V is given by 



dE 
~dt 



, \ ■■ 1 ' i \ -i dV (x) 

TO (X) X + -TO (X) X + 



2 w ax 



^ = °< (25) 



indicating that the total mechanical energy is conserved (unlike equation (10) of Bagchi et al. [16J]) and is the 
constant of motion. Therefore, the PDM-buproducted reaction force Rpdm {%> %) turns out to be a conservative force 
with explicit dependence on both position and velocity. Should we assume that m (x) and m (x) are both positive 
valued functions, then Rpdm (x, x) would act in the direction of the time rate of change of the linear momentum 
(i.e., Rpdm (x,x) = dp/dt = m (x)x 2 /2). The consequences of the very existence of such a force are discussed by 
Mazharimousavi and Mustafa .IT} through the use of the Euler-Langrange's equation (19) and the related examples 
reported therein. 

For the sake of fairness and/or completeness, if one assumes that a PDM particle m (x) = 1/ ' \J\ + £,x 2 is moving 
under the influence of the set of mass- independent forces (as apparently suggested by Bagchi et al. \LM), only then 
Eq.(12) of Bagchi et al. would need a multiplicity of order 1/2 in the first term of y (as a consequence of our proposed 
amendment of Newton's low in (19)). This should not affect their numerical study of the propagation dynamics and 
the sensitivity of the PDM-index £ to enhance phase transition from a limit cycle mode to a chaotic regime and the 
initiation of the complicated nature of bifurcation, etc. 
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